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Calibration of Passive Microwave Polarimeters that 
Use Hybrid Coupler-Based Correlators 

J.R. Piepmeier, Member, IEEE 


Abstract — Four calibration algorithms arc studied for mi- 
crowave polarimeters that use hybrid coupler-based correlators: 
1) conventional two-look of hot and cold sources, 2) three looks 
of hot and cold source combinations, 3) two-look with correlated 
sour ce, and 4) four-look combining methods 2 and 3. The system- 
atic errors are found to depend on the polarimeter component 
parameters and accuracy of calibration noise temperatures. A 
case study radiometer in four different remote sensing scenarios 
was considered in light of these results. Applications for ocean 
surface salinity, ocean surface winds, and soil moisture were 
found to be sensitive to different systematic errors. Finally, a 
standard uncertainty analysis was performed on the four-look 
calibration algorithm, which was found to be most sensitive to 
the correlated calibration source. 

Index Terms — calibration, correlators, microwave radiomefry, 
polarimetry, polarization 


I. INTRODUCTION 

P ASSIVE microwave polarimetry, or polarimetrie radiom- 
etry, is relatively new to Earth remote sensing. While it 
has been quite useful over the last several decades in radio 
astronomy, only recently has microwave polarimetry been 
shown to be of utility to sensing geophysical quantities. Ocean 
surface wind direction remote sensing was indeed enabled 
by polarimetrie radiometers (e.g., [1], [2]). Additionally, the 
use of polarimetry for removing the effects of ionospheric 
Faraday rotation was proposed in [3] and electronic polar- 
ization basis rotation was described in [4]. We define a 
polarimeter, or polarimetrie radiometer, as a radiometer that 
is designed to measure polarization information beyond the 
conventional dual-polarized radiometer. For example, a po- 
larimeter might measure the correlation between the vertically- 
and horizontally-polarized fields, or it might measure ±45°- 
linear polarization. This information is then used to determine 
parameters of the polarization ellipse or the Stokes vector. 
Some examples of Earth science polarimeter architecture can 
be found in [5] with specific implementation examples in [6] 
and [21; for a survey of astronomical polarimetry see [7]. 

In this paper we are specifically studying the calibration 
of a hybrid coupler-based polarimeter of the type shown in 
Fig.l. Our objective is to determine the effects of calibration 
algorithms on the systematic errors of the measurement. The 
polarimeter measures the first three modified Stokes parame- 
ters in brightness temperature T v , T hy and TV. The modified 
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Stokes vector in brightness temperature is defined: 
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where A is the wavelength, k is Boltzmann’s constant, B is the 
detection bandwidth and rj is the intrinsic impedance. The first 
two modified Stokes parameters are the vertical and horizontal 
polarized brightness temperatures measured by conventional 
dual polarization radiometers. The third Stokes parameter is 
measured by correlating the vertical and horizontal voltage 
signals. The hybrid coupler-based polarimeter is based upon 
the principle of correlation through synthesis of the ±45° 
linear polarizations from the vertical and horizontal signals: 

{|2? p | 2 } = I (|£„| 2 ) + \ (|E,| 2 ) +R e(E v E* h ) 

(|£m| 2 ) = \ (|£„| 2 ) + \ (l^l 2 ) - Ke(E v E* h ) 

where the subscripts p and m stand for plus and minus forty- 
five degrees, respectively. 

We do not include the antenna in this analysis. The antenna 
polarization effects are assumed to be compensated by the 
antenna pattern correction process, which is a conventional 
technique used in radio astronomical polarimetry [7]. This 
strategy allows design and calibration methodologies each 
to be separated into two independent parts at the antenna- 
radiometer interface. 

The polarimeter can be described by the following forward 
system equation: 
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where the elements of the left-hand vector are proportional to 
the average voltage at the detector outputs. We can assume 
that the cross-polarization terms in rows one and two of the 
gain matrix are made negligible through the appropriate use of 
isolators and amplifiers. The radiometer contains two amplifier 
and filtering chains that form a conventional dual-polarization 
radiometer used to measured T v and TV In addition, signals 
from the two chains are combined using a 180° -hybrid coupler 
and then detected to measure two additional slant-linear polar- 
ization brightness temperatures T p (for +45° linear or p-pol) 
and T m (for -45° linear or m-pol). The third Stokes parameter 
can be found from their difference: TV = T v - T m . 
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To study and quantify the quality of polarimeter calibration, 
we examine four calibration algorithms using a combination of 
hot, cold, and correlated calibration sources. Each calibration 
methodology produces systematic errors that must be con- 
trolled through design requirements and affect the polarimeter 
calibration uncertainty. The following methodology is used to 
study the four cases. First, a comprehensive forward model 
of the radiometer is developed in Section II. This model is 
assumed to describe how signals actually propagate through 
the polarimetric receiver. Second, for each specific calibration 
case a modified radiometer equation is chosen such that it has 
a reduced number of degrees-of-freedom commensurate with 
the chosen calibration scheme. For example, if only a pair 
of hot and cold calibration looks are provided (two degree- 
of-freedom), then the system can only be modeled with two 
parameters per channel. This modified equation can be inverted 
to arrive at an estimation of TV called TV. Next, calibration 
equations are formed for computing the gains and offsets 
present within the modified radiometer equation. Into these 
equations, expressions for the radiometer outputs are inserted 
based on the actual forward model. Thus, the gains and offsets 
are found in terms of radiometer hardware parameters and 
calibration noise temperatures. Finally, these gains and offsets 
are inserted into the TV equation and the systematic errors are 
found. In Section TV three scenarios (ocean surface salinity, 
ocean surface wind vector, and soil moisture remote sensing) 
using the polarimeter are considered. Given each of these 
three scenarios, quantitative examples of systematic errors are 
given for different calibration algorithms and a quantitative 
uncertainty analysis is completed using the guidelines in [8]. 

II. Polarimeter Forward Model 

In this section, the gain and offset parameters in (3) are 
computed based on the microwave properties of the radiometer 
circuits. Specifically, the parameters can be found in terms 
of the amplitude and phase behavior of the amplifier chains 
and hybrid coupler. This analytically-based model can then be 
used to predict the actual behavior of the radiometer and is 
useful for determining the performance of a chosen calibration 
method, which is done in the remaining sections. 


The amplifier and filter chains have arbitrary transfer func- 
tions #31 (/) and tf 42 (/). The effects of the signal splitters, 
which couple the signals to the hybrid coupler, can be included 
in the transfer functions and do not need to be explicitly noted. 
The hybrid coupler is assumed to have negligible frequency 
dependence and to be well matched over the band of interest, 
which is not unreasonable for operation only over several 
percent of bandwidth. It can be shown that such a hybrid 
coupler, if reciprocal and lossless, has the scattering parameter 
matrix of the form: 

o o s 

•S _ 0 0 vT - sV* -ae** 

s Vl - sV* 0 0 

\yi - s 2 -set* 0 0 / 

( 4 ) 

The reference phase planes are adjusted so that the scatter- 
ing parameter s is real and positive, and lOlogioj-f^ is 
the amplitude imbalance of the hybrid coupler in decibels. 
Amplitude imbalance is plotted versus the scattering parameter 
s (normalized by 2” 1 / 2 ) in Fig.2. Ideally s = 2“ x / 2 and the 
imbalance is 0 dB. The phase term i/> is the phase imbalance 
and ideally is equal to zero. 

The mean voltage out of each detector can be found using 
the above models: 
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where c$ is the detector sensitivity for polarization channel 
V. The signum function in the exponent makes the phase-shift 
factor Hermitian. The gain-bandwidth products G n B, n =1 
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Fig. 2. Hybrid coupler amplitude imbalance in decibels versus the scattering 
parameter s normalized by 2"* 1 / 2 . 


value can be identified: 1) hybrid coupler phase imbalance 
and/or uniform amplifier chain phase imbalance, 2) amplifier 
chain random phase differences, and 3) amplifier chain pass- 
band ripple differences. Below, the parameter a e is computed 
for each of these modes taken independently. First, hybrid 
coupler phase imbalance (ip) and/or constant phase difference 
between amplifier chains (</>) would appear as a cosine factor: 

a e (0, ip) = cos(<p + ip) (10) 

The phase difference between the two amplifier chains 
across the passband, however, is rarely constant. If the phase 
difference is modeled as a uniform random variable with mean 
<p and limits <p±A<p, then a e can be found from the expected 
value of the inner-product integral (9): 

2 He f + °° Hsi (/) H; 2 (/) e^^'df) (11) 

Jo i V 

where the operator is the expected value operator of the 
uniform random variable cp l over the interval ±A(j>. Carrying 
out the integration and factoring out the coefficients identifi- 
able with a e yields: 

a € (<p\ A(p) = cos <p sine A <p (12) 


or 2, are defined as: 
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The gain of channel two can be referenced to channel one 
using a gain imbalance g: 


(?2 = gG\ (7) 

where | log 10 (£) j is the gain imbalance in decibels. The offsets 
resulting from detected receiver noise T n> nx are: 

o v = c v kTi tR xGiB 
on = c v kT 2 ^RxG 2 B 
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The integral in the expressions for the polarimetric channels 
v p and v m is essentially an inner product of H^i and H 42 and 
by the Cauchy-Schwarz theorem, we can write: 

/ -hoo 
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where a e < 1 is a bandpass equalization efficiency. It is 
affected by the differential phase and amplitude variations 
between the amplifier chains and by the phase imbalance of the 
hybrid coupler. A similar expression arises in radio interferom- 
etry and is similarly affected by the channel characteristics [9]. 
It is also effected by decorrelation caused by a group-delay 
difference (so-called fringe washing). Here we assume that 
path lengths are equalized to cancel this effect. 

While a e is determined by the combined effects of phase 
and amplitude variations, three different modes that affect its 


The third mode to examine is passband ripple differences. 
The transfer functions can be parameterized to include sinu- 
soidal amplitude ripple such that 
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where 7 „ is the ripple amplitude and N is the number of 
poles of a representative odd-order Chebyshev filter. The peak- 
to-peak ripple in decibels is 201og 10 | jj^j . (Note that these 
transfer functions still satisfy (6).) By substituting the above 
into the inner-product integral, we find that 
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or if the ripple of the two channels is of equal amplitude j-y), 
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Fig. 3. Dependence of a e ( 7 ) on peak-to-peak amplitude ripple of the 
passband response of the radiometer channels. A 2.5-dB ripple results in up 
to 2% reduction in efficiency. 


Fig. 3 shows the dependence of a e on peak-to-peak amplitude 
ripple in dB. 

Finally, the polarimeter forward model is found by rewriting 
(5) as a matrix equation of the form in (3) and is shown at 
the bottom of the page (17). Throughout the remainder of 
this paper this forward model is used to describe how the 
polarimeter actually operates. It captures nonidealities such as 
differences in gain and phase of the amplifier and filter chains, 
and phase and amplitude imbalance of the hybrid coupler. 
Using this model, we determine the systematic errors of the 
polarimeter under four different calibration algorithms or cases 
in the following section. 

III. Calibration 

Total power radiometers are conventionally calibrated using 
two sources of differing brightness or noise temperatures [10]. 
The polarimeter considered here is essentially four total- 
power radiometers, each sensitive to a different polarization, 
and is tempting to be calibrated as such. The four radiome- 
ters, however, are not entirely decoupled, and the calibration 
scheme can should take this into account. Here, we consider 
four calibration algorithms and examine their effects on the 
systematic errors of the polarimeter. The errors are manifest 
as either or both of polarimetric inefficiency or depolarization 
(see [7]) and cross-polarization coupling depending upon the 


choice of calibration method. Case 1 is the conventional two- 
look method. Case 2 is a three-look method that includes 
varying combinations of hot and cold sources. Case 3 is also 
a three look method but combines Case 1 with a correlated 
calibration source. Finally, Case 4 uses the correlated source 
in conjunction with the three looks of Case 2 for a complete 
characterization of the polarimeter. In all four cases, we 
assume that the two total-power radiometers measuring T v 
and Th can best be calibrated using the conventional two-look 
technique. We, therefore, only consider the calibration of T\j. 


A. Case 1: Conventional hot and cold calibration 

The first case is simply the conventional two-look calibra- 
tion method used with standard total-power radiometers. Each 
channel is calibrated independently using a hot look and cold 
look by viewing two calibrators with differing brightness or 
noise temperatures. Because only two calibration states are 
used, only two calibration parameters can be found for each 
channel. The following is used for the Case-1 radiometer 
model for measuring T p and T m : 



As is seen above, only two parameters are needed for each 
channel: G p and o p for the +45° polarization (p-pol) and G m 
and o m for the -45° polarization (m-pol). Solving for Tu 
yields the estimate Tjy: 


Tu = 
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The two calibrators are implemented as a matched load for 
the cold noise temperature T c and an additional coupled noise 
diode for the hot noise temperature Th as seen in Fig. 1. For 
the p-pol channel, the calibration equation is 
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\ v P,nJ V» V \Op) 


or solving for G p and o p 


( 20 ) 
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The noise temperatures Tc and Th are the nominal noise 
temperatures (including errors) of the calibration sources. 

Now, if the actual observed voltages as calculated by the 
complete forward model (17) are substituted into the gain and 
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offset equations, the estimate Ty can be found in terms of the 
input Stokes vector. First, we find that 

v p , c = c p kB(s 2 G 1 + (l-s 2 )G 2 )T^ _ 

v Pi H — c p kB {s 2 Q x + (l - s 2 ) G 2 ) V H 

where the primed temperatures indicate the true noise tem- 
peratures of the calibration sources). Substituting this into the 
gain and offset expressions: 

G p = c p kB(s 2 G 1 + (l-s 2 )G 2 )%^ 

o P = CpkB^GiTi'jix + (?2 (l — s 2 ) T 2 ,rx (23) 

+ (s 2 G 1 + (l-s 2 ) g 2 ) 

and similarly for the m-pol channel: 

G m = C m kB ((1 - S 2 ) Gl + 3 2 G 2 ) ^ 

— S 2 ) G\.T\.RX + S 2 G 2 T 2 'RX 

(24) 

Going the next step, substitute these into Ty: 




y(g(l+g)avT^ r g -r c 
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Here the inadequacies of conventional two-look calibration are 
seen quite clearly. First, there is a gain term affecting Ty, part 
of which is directly proportional to the bandpass equalization 
efficiency a e . As will be seen in Case 3, this error can be 
removed by using a correlated calibration source. Amplitude 
imbalance in both the receiver channels and the hybrid coupler 
also affect the gain of Tu or the polarimetric efficiency, as it 
is called in the astronomical polarimetry literature. 

The other error is cross polarization contamination. Ampli- 
tude imbalance in the hybrid coupler actually causes mixing 
of the second (not modified) Stokes parameter ( Tq = T v -TO 
into Ty. This is easily shown by taking (25) and setting all 
parameters to their ideal values (e.g., g = 1) except for s, 
which yields: 

fu(s) = (2s 2 - 1) (T v - T h ) + 2s\/l - s 2 Tu (26) 

Fig. 4 shows the amount of Tq mixing into Ty as a function 
of hybrid coupler amplitude imbalance. There is about 11.5% 
mixing per decibel of imbalance. The next calibration method 
Case 2 uses independent hot and cold looks to remove cross- 
polarization mixing. 


B. Case 2: Independently applied hot and cold looks 



Fig. 4. Mixing ratio of Tq into Tu due to hybrid coupler imbalance for 
Case 1 (conventional two-look calibration). 


and cold looks, the third state is one channel viewing hot with 
the other viewing cold. This third degree of freedom allows 
the radiometer model to be slightly more complex than (18). 
The case-2 radiometer model is: 

( V p\ = - (G pv \j GpyGph \ f j? \ ( °p\ 

\VmJ 2 \(? mv Gjjih —y/G pv GphJ l rp^ J \OmJ 

(27) 

The differentiation of the gain terms explicitly includes the 
contributions of the vertical and horizontal polarizations to the 
p-pol and m-pol channels. Solving for the radiometer equation 
to find the estimate Ty: 


Tu = ■ g ~ ( Vp ~° p ~ fhGph ~ fpGp ) 

v'mh^mv i ^ph^pv ' ' 

VGmhG mv 

G m hG mv + GphG vv 


(fim °m ThGjnfi T v G m -^j 


(28) 


where the estimates T v and T h are taken from their respective 
total power channels. To take advantage of this type of 
calibration method, the calibration topology of the radiometer 
must allow for the independent switching of the hot and cold 
sources. This is readily done by literally implementing the 
configuration shown in Fig.l. One such way is to use a switch, 
a matched load, and a coupled noise diode. On the other hand, 
if the calibration was thru the feedhom, then some sort of 
externally polarized source would be needed and the explicit 
separation of the antenna from the receiver would be violated, 
further complicating the system. 

The gains and offsets are found as in Case 1: 


In situations where the cross-polarization mixing effect / u P) c \ (Tc Tc l\ ( G pv \ 

dominates the overall uncertainty (as in (26)), a third cali- v P: h I = l Ty Th 1 G pm J 

bratiori state can be added to Case 1. In addition to the hot \Vp,CH J \Tc T H 1 / \ o p J 


( 29 ) 
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which when solved yields: 


r _ V P>H ~ V P;CH 

PV T H -Tc 

r - V ?' CH ~ v p> c 
ph T h -T c 


(30) 

(31) 


Inverting this model results in the following estimate of TV: 
Grn (u p Op) G p (f m — O m ) 


Tu = 


(35) 


and Op is the same as for Case 1. The m-pol coefficients are 
can be found in a similar fashion. 

Following the same steps performed in the previous section 
results in the following expression for TV: 


G m G pU GpG-mu 

The gains and offsets are found as before, using p-pol as an 
example: 


(36) 




The third calibration state primarily compensates for cross 
polarization mixing. In essence, any residual offset is due to 
errors in estimating T v and T h . For example, if we simplify 
(32) by assuming that the amplifier chains are gain balanced 
(g = 1), the detectors have the same sensitivity (c m = c p ), 
and that the errors in estimated noise temperatures of the 
calibration looks are negligible, then 

fu = CC'TU + (Tv -T v -T h + T h ) (33) 

where T v and Th are the true brightness temperatures and 
T v and T h are the quantities estimated using the radiometer. 
Note that the residual offset is of similar form to the cross- 
polarization coupling in Case 1; however, the magnitude is not 
proportional to Tq , but rather the error in the estimate of Tq. 

As in Case 1, the bandpass equalization efficiency and errors 
in the hot and cold noise temperatures act on the polarimetric 
efficiency. In fact, their effect is identical to Case 1. Unlike 
Case 1, however, the third calibration state has removed the 
effect of gain imbalance g and coupler amplitude imbalance s 
on the polarimetric efficiency The remaining effects only deal 
in phase. So for a given set of hardware, the third calibration 
state slightly improves the polarimetric efficiency in addition 
to nearly complete compensation of cross polarization mixing. 


C. Case 3: Calibration with correlated or polarized source 


On the other hand, in situations where polarimetric effi- 
ciency is the dominating consideration, as opposed to cross 
polarization coupling, the third calibration state can be re- 
placed a with correlated source. Thus, the gain of the third 
Stokes parameter channel can be explicitly included in the 
case-3 radiometer model: 




The value of the correlated noise temperature Tcn depends 
upon implementation. For example, if the output of a noise 
diode is split and then coupled into the two channels, Tcn 
is the diode noise temperature reduced by the coupling factor 
of the directional coupler. The sign of T C n is determined 
by the relative phase of the coupled signals. When the noise 
diode is turned on, it adds power to the total power channel 
and generates a correlated signal mimicking a third Stokes 
parameter input. Doing so allows G p u to be found: 


G p u = 


Vp,CN ~ Vp,C 

Tcn 



(37) 


where G p is given by 21. In addition to that for G p , the result 
for o p is identical to Case 1. The m-pol results are similar. 

After substituting the calibration coefficients into the ra- 
diometer equation, the following is found for TV: 


Tu = 


Tcn T V? 2s 2 — 1 T C n -i 
Tcn (1 +p) sV 1 — $ 2 Tq N 


(T v - T h ) (38) 


Now the calibration of TV is directly dependent upon the ac- 
curacy with which the correlated noise temperature is known. 
The polarimetric efficiency is compensated by the correlated 
noise source and any error in Tcn is transferred directly to 
TV- Additionally, like Case 1, Tq is mixed directly into TV by 
coupler imbalance. The mixing ratio is within 0. 1 percentage 
points of the case-1 result shown in Fig.4. 


D. Case 4: Full calibration 

In the final case, the calibration schemes of the previous 
cases are combined to fully calibrate the radiometer. In other- 
words, four calibration states are used: conventional hot and 
cold (from Case 1), mixed hot and cold (from Case 2), and 
correlated noise (from Case 3). With four calibrators come four 
degrees-of-freedom, allowing the polarimeter to be described 
by the following forward system equation: 

(Z) = (oZ al oZ) (|) + (Z) (39 > 

Note, the above is exactly the bottom two rows of the complete 
model in (3). The third Stokes parameter estimate using this 
system is: 

Tu — , y-f2 ( Vm ~ ° m ~ T v G mv — ThG m h\ 

U mU + ij P U V ' 

+ Qi q2 ( v p ~ °p " T v G pv — ThGphj (40) 
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Three gains and one offset can be found for each channel. 
Again using the p-pol as an example, the calibration equation 
is: 



T C 0 1\ 


G pv \ 

T h 0 1 


Gpm 

Th 0 1 


G p u 

Tc + \Tcn ±Tcn 1/ 


\°p J 


(41) 


by careful amplitude and phase balancing of the radiometer. If 
we assume, however, that all noise and brightness temperatures 
are known without error, then (42) degenerates to the error free 
expression Tu = T\j. Thus, there is a trade-off between the 
accuracy of the calibration noise temperature values and the 
degree to which the radiometer components are balance. Of 
course it is not possible to know the calibration and brightness 
temperatures with complete accuracy, so one must consider 
their effect. 

Consider the simplification where the radiometer is perfectly 
balanced but the calibration noise temperatures and T v and T h 
estimates are only known with error: 


Not surprisingly, the expressions for the first two gain coeffi- 
cients G pv and G mv are the same as for Case 2, G p u is the 
same as Case 3, and the offset is the same as Case 1. Similar 
results occur for the m-pol channel. This relation is indicative 
of the effects of the different calibration sources on the overall 
instrument calibration. The hot and cold looks help identify the 
offset, adding the hot-cold combination look differentiates the 
cross-polarization gain coefficients; and finally the correlated 
source, the third Stokes parameter gain coefficient. 

The resulting estimate for Tu is quite lengthy. It is, however, 
conceptually quite simple: the accuracy of the estimate is 
affected by how well the calibration noise temperatures are 
known and how well the radiometer channels are balanced in 
gain and phase. The full expression is: 
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As in the previous cases, we can identify polarimetric effi- 
ciency and cross-polarization coupling errors in the expression 
for Tj. In the error terms there is an interaction between the 
radiometer parameters. On one hand, systematic errors in the 
measurement are caused by errors in calibration noise temper- 
ature values, but the effect of these errors can be minimized 


T V = 


Tcn 


T C n H 




where 




A T l — ' 

ThT’c - T C T' H 

Th-Tc 


(43) 


H= 1- 


T' h - T' c T cn 
Th - T c Tfr N 


There are basically two main error modes in the above. First, 
any error in the nominal correlated noise temperature Tcn is 
transferred directly to Tu. This dependence is not unexpected 
because the correlated noise source was added to compensate 
for the polarimetric efficiency. Second, errors in the hot and 
cold noise temperatures and/or T v and T h cause offsets by 
altering the contribution of T v and T h to the estimate. At this 
point, it is more instructive to examine the errors quantitatively 
by considering a case study radiometer. 


IV. Applications 

In the previous section we have shown how the combination 
of imperfect hardware components and selection of calibration 
method causes systematic errors. In this section the results 
are applied to several remote sensing problems quantitatively 
demonstrating the effect of choice of calibration algorithm. 
Four scenarios are considered: ocean surface salinity (OSS), 
ocean surface winds (OSW), and soil moisture (SM) with 
and without antenna polarization basis rotation. In both OSS 
and SM remote sensing the third Stokes is used for Faraday 
rotation correction [3], [11]. In OSW, the third Stokes param- 
eter is used in wind direction retrieval, e.g. [1], For the SM 
scenario, we add an additional case SM(b) that considers a 
45° polarization basis rotation within the antenna system. Such 
is relevant to a scanning radiometer antenna that has a fixed 
feed and a rotating reflector where the polarization basis is 
continuously rotated and then electronically corrected in post 
processing [4]. 

We concentrate on examining the systematic errors and cal- 
ibration uncertainty arising from systematic effects. We define 
this component of uncertainty based upon practices of the U.S. 
National Institute of Standards and Technology (NIST) [12] 
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TABLE I 

PROPERTIES OF HYPOTHETICAL CASE STUDY RADIOMETER 


Component parameters 


channel gain imbalance 

2 dB 

channel phase imbalance 

10° 

channel phase variation 

20° 

channel amplitude ripple 

2 dB 

coupler amplitude imbalance 

0.2 dB 

coupler phase imbalance 

5° 


Resulting model parameters 


s 

0.700 (from Fig, 2) 

9 

1.585 

«e(7) 

0.987 (from Fig.3) 


0.946 from (12) 

Qc « a e (7 + 

0.934 


Calibration parameters 


T C 

250 K 

Th 

350 K 

Ton 

50 K 


and International Standards Organization (ISO) [8]. A sys- 
tematic effect is something that can produce a systematic 
error (or bias) in the measurement process. For example, 
uncertainty in an estimated calibration target emissivity would 
contribute to the uncertainty arising from systematic effects 
for a standard microwave radiometer. On the other hand, 
the random noise in a radiometer measurement due to finite 
integration time (i.e., NEAT) does not fall under this category 
of uncertainty. Another way to consider this uncertainty is 
that which cannot be reduced by increasing the number of 
measurements. In our situation, the systematic effects and 
errors and their associated calibration uncertainly are driven by 
gain and phase balance behavior of the amplifier chains and the 
hybrid coupler and choice of calibration method. Performing 
multiple measurements or increasing integration time is not 
going to reduce the uncertainty caused by these effects. 

The properties of a hypothetical case-study radiometer of 
the form in Fig. 1 are listed in Table I. These values where 
chosen to represent some typical performance specifications 
of commercial off-the-shelf (COTS) components. Of course 
these values could differ depending upon the radiometer 
center frequency, bandwidth, and amount of effort invested 
in design and construction. Nonetheless, they are useful for 
demonstration purposes. 

For the first application, the estimated Tu is computed 
for each of the four different measurement scenarios using 
the first three calibration algorithms. We assume that the 
calibration noise temperatures and T v and T h are perfectly 
known. This allows us to demonstrate numerically the effect 
of hardware nonidealities and different calibration methods on 
the measurement result. (Calibration Case 4 was not included 
because there are no systematic calibration errors under our 
assumptions.) The results are presented in Table II. A general 
description of the scene, the incidence angle, observation band, 
and brightness temperatures are listed in the top rows of 
the table. The brightness temperatures were determined using 
published data [10, Figs. 18.2, 18.19, 18.20 and 19.12]. For 
each calibration case and remote sensing scenario the values 


TABLE H 

Third stokes Parameter Estimates and Case Study results 


Scenario 

OSS 

OSW 

SM (a) 

SM (b) 


Scene 

Ocean 

Ocean 

Land 

Land 


Circumstance 

12° IFR 

Winds 

7° IFR 

45° PBR 


Incidence angle 

33° 

53° 

40° 

40° 


Frequency band 

L 

K 

L 

L 


T v 

105 

180 

215 . 

198 

K 

t k 

80 

120 

170 

OO 

oo 

K 

Tu 

10 

0.5 

10 

-45 

K 

Tu 

8.63 

-0.68 

8.25 

-41.1 

K 

error 

-1.37 

-1.18 

-1.75 

3.87 

K 

Case 1 , 

(CC HH) rC ' CITOr 

-14% 

-236% 

-18% 

-8.6% 


gain (m) 

0.91 

0.91 

0.91 

0.91 


offset (6) 

-0.47 

-1.13 

-0.85 

-0.19 

K 

Tu 

9.34 

0.47 

9.34 

-42.0 


Case 2 error 

-0.66 

-0.03 

-0.66 

2.96 

K 

(CC,HH, rel. error 

-6.6% 

-6.6% 

-6.6% 

-6.6% 


gain (ttv) 

0.93 

0.93 

0.93 

0.93 


offset (b) 

0 

0 

0 

0 

K 

T V 

9.48 

-0.75 

9.07 

-45.2 

K 

Case 3 error 

-0.52 

! -1.25 

[ -0.93 

-0.21 

K 

(CC,HH, rel. error 

-5.2% 

-249% 

t -9% 

-0.5% 


CN) gain (m) 

1 

1 

1 

1 


offset ( b ) 

-0.52 

-1.25 

-0.93 

-0.21 

K 


IFR=ionospheric Faraday rotation 
PBR=polarization basis rotation 


for Tu and associated absolute and relative errors are listed. 
While under our assumptions expressions for Tu (25), (32), 
(38) and (42) can be summarized by 

Tjj — mTu + b (44) 

where m is the Tu gain and b the Tu offset. Values for m 
and b are given in Table n as well. These values are useful for 
determining the effects of the different calibration algorithms 
upon the errors. 

For OSS, the errors are all >0.5 K, which is higher than 
desired for the Faraday rotation correction application. An 
error budget of ~0.1 K is stated in [3]. Thus, the case-4 
calibration algorithm is most likely needed. Clearly cross- 
polarization errors are detrimental for OSW. In OSW the 
polarization contrast Tq is relatively large (~60 K) and cou- 
ples into Tu via imperfect hybrid coupler amplitude balance. 
Because Tu signatures of OSW are quite small «1 K), any 
bit of Tq coupling into Tu is detrimental. The error is slightly 
worse for Case 3 because the Tq leakage is emphasized by 
the polarimetric efficiency correction that occurs when the 
correlated noise source is included in the calibration. 

The errors for SM(a) are similar to the case for OSS. The 
required brightness temperature uncertainty for SM, however, 
is at least an order-of-magnitude less stringent than for OSS, 
which will relax the Tu calibration requirement. On the other 
hand, the absolute errors for SM(b) in Cases 1 and 2 are much 
larger than the other scenarios. The larger value of Tu coupled 
with a Tu gain of ~0.9 causes ~ 3-K error in Tu, Case 
3, however, readily demonstrates the utility of the correlated 
noise source; the gain is corrected reducing the absolute error 
to a residual 0.2 K, which is from cross-polarization coupling. 



PREPARED FOR IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING, PRINTED ON MARCH 31, 2003 


9 


TABLE m 

Uncertainty Analysis Results - Case 4, OSS Scenario 


Sensitivies (K/K) 

| Uncertainties (K) 

Products (K) 

dT v /dT H 

-0.0216 

ut h 

0.5 

0.0108 

dTu/dTc 

0.0315 

U T C 

0.5 

0.0157 

dTu /dTcN 

0.2010 

Ut CN 

0.5 

0.1005 

OTu/dTv 

0.0169 

u f„ 

0.5 

0.0085 

dTrj/dn 

-0.0268 

. "4- 

0.5 

0.0134 


Total combined uncertainty by R.S.S. (K) 0.1035 


Then sensitivity is driven by the ratio of the value of T v to 
the value of correlated noise temperature. Having a larger 
value of correlated noise temperature would improve the 
calibration because the errors are interpolation errors when 
\Tu\ lies between 0 and Tcn • Of course, this only holds if 
increasing the value of Tqn does not increase its uncertainty 
by a defeating amount. The principle is similar to two-look 
calibration of a conventional total-power radiometer, whereby 
an optimal hot load temperature can be found by trading the 
uncertainty and sensitivity [13]. 


A correlated noise source is probably needed when an antenna 
with a rotating polarization basis is used in the radiometer. 
Whether or not correction of Tq coupling using independent 
hot and cold looks is needed depends upon the expected 
amounts of ionospheric Faraday rotation and the overall error 
budget. Further investigation is needed before a conclusion 
can be made as to the optimal calibration algorithm (optimal 
in the sense of the least complex that meets requirements) for 
SM. 

Because Case 4 has no systematic errors when the calibra- 
tion noise temperatures are perfectly known, the simple nu- 
merical examples from above are not instructive. To generate 
errors using Case 4, we must assume that there is uncertainty 
in the calibration noise temperatures. This is typically done 
using the standard combined uncertainty expression: 


u *v = 


dT v 


dTu) 2 

— 1 *4 C 


8T b ) Uth + \dT C/ 

+ (H;) “k- 


+ 


dTu 
df v J 


4.+ 


dtu\ 

dt h ) ^ 




(45) 


where m (i ~ T H ,T C , etc.) are the standard uncertainties of 
the noise and brightness temperatures. Using the numerical 
values for the scene from the OSS scenario, the above terms 
are evaluated and presented in Table m. We use the same 
0.5-K value for all the calibration temperature uncertainties, 
which could arguably be too low or too high for any specific 
uncertainty. Nonetheless, using the same value is important so 
that the different contributions are not unduly weighted for the 
purpose of the example. 

Of significant note in Table m is that the total uncertainty 
is dominated by the sensitivity to the correlated noise tem- 
perature. Its contribution is approximately ten times greater 
than those of the other sources. The uncertainty could thus be 
approximated: 




(46) 


where by examining (42) the sensitivity is approximately 


dTu Tu 
dT CN * T' cn 


(47) 


V. Conclusions 

We have proposed four calibration algorithms for microwave 
polarimeters that use hybrid coupler-based correlators. Be- 
cause amplitude and bhase imbalance in the radiometer com- 
ponents cause systematic errors in the third Stokes parameter 
measurement, careful consideration of calibration algorithms 
is necessary for reducing these errors. The four algorithms 
considered were (1) conventional two-look of hot and cold 
sources, (2) three looks of hot and cold source combinations, 
(3) two-look with correlated source, and (4) four-look com- 
bining methods 2 and 3. Each of the four methods produce 
differing amounts of systematic error. When calibrating with 
method 1, the polarimeter has both gain and offset errors. 
Method 2 removes the offset error, while method 3 the gain 
error. Finally, method 4 compensates for all errors. The results 
obtained for the calibration methods where then applied to a 
case study radiometer in four different remote sensing sce- 
narios. Applications for ocean surface salinity, ocean surface 
winds, and soil moisture with and without antenna feed polar- 
ization basis rotation were considered. This numerical study 
illustrated that the required calibration technique depends upon 
the application. Offset errors were more detrimental than gain 
errors for ocean surface winds remote sensing. Ocean surface 
salinity requires full calibration and soil moisture requires 
at least the correlated calibration source when used with a 
polarization basis rotation antenna system. Because the four- 
look calibration method compensates for all systematic errors, 
a standard uncertainty analysis was performed to elucidate the 
dependencies on calibration source uncertainties. The largest 
sensitivity is to the correlated noise temperature. In particular, 
the value of the noise temperature and its uncertainty can be 
traded to achieve a minimum total uncertainty. 

References 

[\] J. R. Piepmeier and A. J. Gasiewslri, “High-resolution passive polari- 
metric microwave mapping of ocean surface wind vector fields,” IEEE 
Trans. Geosci. Remote Sensing , vol. 39, no. 3, pp. 606-622, 2001. 

[2] S. Yueh, W. Wilson, F. Li, S. Nghiem, and W. Ricketts, “Polarimetric 
measurements of sea surface brightness temperatures using an aircraft 
K-band radiometer,” IEEE Trans . Geosci . Remote Sensing , vol. 33, no. 

1, pp. 85-92, 1995. 

[3] S. H. Yueh, “Estimates of Faraday rotation with passive microwave 
polarimetry for microwave remote sensing of Earth surfaces,” IEEE 
Trans . Geosci. Remote Sensing, vol. 38, no. 5, pp. 2434-24388, 2000. 

[41 A. J. Gasiewski and D. B. Kunkee, “Calibration and applications of 
polarization-correlating radiometers,” IEEE Trans . Microwave Theory- 
Tech vol. 41, no. 5, pp. 767-773, 1993. 

[5] N. Skou, B. Laursen, and S. Sobjaerg, “Polarimetric radiometer 
configurations: potential accuracy and sensitivity,” IEEE Trans. Geosci 
Remote Sensing, vol, 37, no. 5, pp. 2165-2171, 1999. 



PREPARED FOR IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING. PRINTED ON MARCH 31, 2003 


10 


[6] J. R. Piepmeier and A. J . Gasiewski, “Digital correlation polarimetry: 
Analysis and demonstration,” IEEE Trans . Geosci. Remote Sensing, vol. 
39, no. 11, pp. 2392-2410, 2001. 

[7] Jaap Tinbergen, Astronomical Polarimetry , Cambridge University Press, 
New York, 1996. 

[8] International Standards Organization (ISO), Guide to the Expression of 
Uncertainty in Measurement , Switzerland, 1995. 

[9] A. Richard Thompson, James M. Morgan, and George W. Swenson, Jr., 
Interferometry and Synthesis in Radio Astronomy, John Wiley 8c Sons, 
Inc., New York, 2 nd edition, 2001. 

[10] F. T. Ulaby, R. K. Moore, and A. K. Fung, Microwave Remote Sensing: 
Active and Passive , Addison- Wesley Pub. Co., Reading, Mass., 1986. 

[11] D.M. Le Vine and S. Abraham, “The effect of the ionosphere on remote 
sensing of sea surface salinity from space: absorption and emission at L 
band,” IEEE Trans . Geosci Remote Sensing, vol. 40, no. 4, pp. 771-782, 
2002. 

[12] Barry N. Taylor and Chris E. Kuyatt, “Guidelines for evaluating 
and expressing the uncertainty of NIST measurement results,” Tech. 
Rep. 1297, National Institute of Standards and Technology, Gaithers- 
burg, MD, September 1994, Available on the World Wide Web at 
http://physics.nist.gov/Pubs/guidelines/contents.htmL 

[13] P. Racette and R. Lang, “Comparative analysis of radiometer systems 
using non- stationary processes,” in Proc. Int. Geosci . and Remote Sens. 
Symp. (IGARSS), Toronto, Canada, 2002, pp. 565-567. 


Jeffrey R. Piepmeier (S’89-M’99) received the 
Ph.D. degree in electrical engineering from Geor- 
gia Institute of Technology in 1999. Previously, he 
received the M.S. degree in electrical engineering 
from Georgia Tech in 1994 and the B.S. degree in 
engineering from LeTourneau University, Longview, 
Texas in 1993. From 1993 to 1994 he was a Schak- 
leford Fellow with the Georgia Tech Research Insti- 
tute. Piepmeier was third place winner in the 1998 
IGARSS student prize paper competition. In 1999, 
he joined the Microwave Instrument Technology 
Branch at NAS As Goddard Space Flight Center, Greenbelt, MD. Dr. Piepmeier 
was recipient of an Excellence in Federal Career Gold Award (Rookie-of- 
the-year) in 2000 and is a 2002 NASA Earth Science New Investigator. 
In 2000 he served as conference chair of the First Microwave Radiometer 
Calibration Workshop (MicroCal2000). He has been co-chairman of the 
Instrumentation and Future Technologies technical committee of the IEEE 
Geoscience and Remote Sensing Society since 2000. He is a member of 
IEEE, URSI (Commission F), and the American Geophysical Union. His 
technical interests include passive remote sensing and technology development 
for microwave remote sensing. 



